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Abstract 

We consider the planar MHV amplitude in Af = 4 supersymmetric Yang-Mills theory for 
2 — !• 4 particle scattering at two and three loops in the Regge kinematics. We perform an analytic 
continuation of two-loop result for the remainder function found by Goncharov, Spradlin, Vergu 
and Volovich to the physical region, where the remainder function does not vanish in the Regge 
limit. After the continuation both the leading and the subleading in the logarithm of the energy 
terms are extracted and analyzed. Using this result we calculate the next-to-leading corrections 
to the impact factors required in the BFKL approach. The BFKL technique was used to find 
the leading imaginary and real parts of the remainder function at three loops. 
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1 Introduction 



Recently we have witnessed revolutionary developments in studying scattering amplitudes in 
supersymmetric theories. The present progress is traced back to the work of Parke and Tay- 
lor PP, who showed that the tree-level gluon scattering amplitude can be written in a very 
compact form for some particular helicities of the external particles, namely the maximally he- 
licity violating (MHV) amplitude. The simplicity of the Parke- Taylor tree-level formula raised 
a hope that the quantum corrections could be also compactly encoded in the MHV gluon am- 
plitudes. A great effort in this direction led to formulation by Anastasiou, Bern, Dixon and 
Kosower (ABDK) [2] and then by Bern, Dixon and Smirnov (BDS) [3] ansatz for multi-loop 
planar gluon MHV amplitude in JV = 4 super Yang-Mills theory. 

The BDS formula was tested by Alday and Maldacena [1] from strong coupling side using 
conjectured AdS/CFT correspondence in the limit of large number of external legs. They 
argued that the BDS ansatz is probably to be violated starting at six external gluons. This 
violation was established by Bartels, Sabio Vera and one of the authors (BLS) [3] analyzing the 
analytic structure of the BDS amplitude. It was shown that the BDS ansatz for six-particle 
amplitude at two loops is not compatible with the Steinmann relations [6], that impose the 
absence of the simultaneous singularities in the overlapping channels. They also showed [7] that 
the BDS violating piece originates from the so-called Mandelstam cuts, which are the moving 
Regge singularities in the complex angular momenta plane. The BDS violating term in the 
multi-Regge kinematics was explicitly calculated [7] with logarithmic accuracy in the physical 
region, where it gives a non- vanishing and pure imaginary contribution. We call this region the 
Mandelstam region (channel). The BDS violating term for the six-point planar MHV amplitude 
was found using the BFKL approach [5] and for an arbitrary number of external gluons it 
contains contributions of Mandelstam cuts constructed from an arbitrary number of reggeized 
gluons for the Bartels-Kwiecinski-Praszalowicz (BKP) state [9] [TU] with the local Hamiltonian 
of an integrable open Heisenberg spin chain [11] . Other interesting limits of MHV amplitudes 
were studied in the Regge kinematics by Brower, Nastase, Schnitzer and Chung-I Tan [12]. 

Drummond, Henn, Korchemsky and Sokatchev [13] analyzed the conformal properties of 
polygon Wilson loops in M = 4 SYM and showed that anomalous conformal Ward identities 
uniquely fix the form of the all-loop 4- and 5-point amplitudes, so that any relative correction 
to the BDS ansatz starting at six external particles should be a function of conformal invariants 
(cross ratios of dual coordinates). The relative correction to the BDS formula was named the 
remainder function for an amplitude with L loops and n external legs, and the first non- 

(2) 

trivial remainder function is R 6 . 

It was sugg ested QH EH QSH2 US US that R<h L) can be obtained from the expectation value 
of the light-like polygonal Wilson loops. Del Duca, Duhr and Smirnov [2UI I2T] expressed R 6 ~ 
in terms of generalized polylogarithms, which was greatly simplified by Goncharov, Spradlin, 
Vergu and Volovich (GSVV) [22], and then written only in terms of Life functions with arguments 
depending on three dual conformal cross ratios. 

The two-loop remainder function i?g for the scattering of eight gluons was calculated by 
Del Duca, Duhr and Smirnov [23] and its diagrammatic structure was analyzed by Alday [24] , 

(2) 

The form of i?g is remarkably simple and it is constructed only of a product of some logarithms 
plus a constant term. 

Earlier we performed [25] an analytic continuation of the GSVV formula to a physical region 
considered in refs. [5] [7]. The continuation showed a full agreement between the BLS formula 
and the Wilson loop calculations at the leading logarithmic level and allowed to extract the 
terms subleading in the logarithm of the energy. Numerically, an agreement between the two 
approaches was demonstrated by Schabinger [26] . The analytic continuation to the mentioned 
above physical region in the regime of the strong coupling constant was performed by Bartels, 
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Kotanski and Schomerus [27] . They found the leading singularity, which governs the high energy 
behavior of the scattering amplitude, so-called reggeon intercept. At weak coupling constant 
the corresponding intercept is determined by the BFKL equation [TJ. 

In this study we present some details of the analytic continuation performed by the authors 
in ref. [25] ■ Based on the obtained result we calculate the next-to-leading (NLO) impact factors 
for the color octet states in the BFKL approach. In the BFKL technique we also find the three- 
loop contribution to the remainder function of planar six-point MHV amplitude in the leading 
logarithmic approximation (LLA) as well as the real part of the subleading corrections in the 
next-to-leading logarithmic approximation. 

2 BFKL approach 

In this section we briefly outline the results of the BFKL approach to the planar MHV 
amplitudes in M = 4 SYM. 

We consider the six-point MHV amplitude for production of two gluons with momenta k\ 
and hi in small angle scattering of the particles with momenta pa and as depicted in Fig. [TJ 




DiS % M 2->4 D 'S\ M 2 _ >4 

(a) m = e l0 |iti| (s, sa,si,S3, soi2>si23 > 0) (b) in = e~ l27T \m\ (s, s 2 > 0; si, s 3 , floia, S123 < 0) 



Figure 1: The BDS violating contribution appears in the region s, S2 > 0; si, S3 < 0. 
All energy invariants are shown in Fig. Q] and are related to the dual conformal cross ratios 



by 



s s 2 S 1 t 3 S3 ti 

Ui = , u 2 = — , u 3 = — . (1) 

S012 S123 Sq12 E2 S123 12 



The multi-Regge kinematics is equivalent to having s ^> soi2,si23 3> 81,82,83 3> ti,t2,ta, 
which in the terms of the cross ratios reads 



1- «i -►+(), u 2 ^+0, w 3 ^+0, —=—siO{l), — 2-~0(l). (2) 

1 — Ui 1 — Ul 

In this kinematics the remainder function of the MHV amplitude goes to zero in direct 
channel in Fig. QJi, while in the Mandelstam channel Fig. [TJd grows with In S2 and becomes pure 
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imaginary. In the Mandelstam channel the gluon momenta k\ and &2 are flipped and the cross 
ratio u\ possesses a phase 



u lb =u la e-^, (3) 

leading to necessity of an analytic continuation of Rq. It was demonstrated by Bartels, Sabio 
Vera and one of the authors [7] that the BDS violating piece comes from the Mandelstam 
cut state propagating in the crossing channel between the produced particles k\ and fe, and 
denoted by the dark box in Fig. [1] In J\f = 4 SYM (as well in QCD) for a large number 
of colors this state is described by the color octet BFKL evolution equation [8]. The BFKL 
equation can be formulated as the Schrodinger equation with a Hamiltonian equivalent to that 
of a completely integrable open Heisenberg spin chain model [H] , which made it possible to solve 
it analytically [7] . In the direct channel of the multi-Regge kinematics (see Fig. QJ,) given by 
Eq. [5] the remainder function vanishes due to the Mandelstam cancellation of cut contributions 
as was shown in ref. [S]. 

The BDS violating piece in the Mandelstam channel is given by [7] 

M 2 ^ 4 = M 2 B f 4 s (l + zA 2 ^ 4 ), (4) 

where is the BDS amplitude [3] and the correction A2_>.4 was calculated in all orders with 

a leading logarithmic accuracy using the solution to the octet BFKL equation. The all-orders 
LLA expression for A2_j.4 reads 

— !„£ <-)"/:^? (tr (&p (^-h < 5 > 



/ (_i)« r dv 2 ffiV" 7 ( ( i _ Ul) -Kn) _ a 



n— — oo 



Here &i , &2 are transverse components of produced gluon momenta, qi , qi , 93 are the momenta 
of reggeons in the corresponding crossing channels and 

u{v, n) = -aE^ n . (6) 

The perturbation theory parameter a and the eigenvalue of the color octet BFKL E v n are given 

by 

- - 9%K (Ane-^Y ^ (7) 



8tt 2 v ' 2ir 



and 



Eu, n = + * (1 + + y) + ^ ( x - w + t) " 2 ^ (1) ' (8) 

where tf)(z) = r'(z)/r(z), 7 is the Euler constant 7 = —ip(l) and the dimensional regularization 
parameter e is defined by d = 4 — 2e. 

The second line of Eq. [5] follows from the fact that in the Regge kinematics the energy 
invariant S2 is related to the cross ratio u\ by 

l- Ul ^±±hl= s J±. (9) 

This way we set Sq = (ki + k2) 2 to be an energy scale, which becomes relevant only beyond 
leading logarithmic approximation. The choice of the energy scale sq is natural in the Regge 
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kinematics because it reflects the smallness of the transverse components with respect to the 
longitudinal components of the particle momenta. The expression in Eq. [5] is a function of the 
dual conformal cross ratios Ui as discussed in section |U The BDS violating piece at two loops 
found in ref. [7] can be written in terms of the reduced cross ratios as 

using Eq. [§] and the fact that the reduced cross ratios 

"2 = , «3 = Z (11 

1 — u\ 1 — ui 

in the multi-Regge kinematics are given by 

^ |k 2 | 2 | qi | 2 ^ / 19 x 

" 2 " |k 2 + k 1 |W ^ 3 — |k 2 + k i p|q 2 p " M 

Surprisingly, the expression in Eq. [TU]can be obtained from the BDS formula using only gen- 
eral analytic properties of the scattering amplitudes and the factorization hypothesis (proposed 
by Alday and Maldacena [14]) as it was shown by one of the authors [28 . In this technique 
it is enough to know the form of the BDS amplitude at one loop to obtain the leading loga- 
rithmic imaginary term at two loops. Unfortunately, for the three loops the knowledge of only 
the BDS formula is not enough and some extra information is to be included in the analysis. 

(2) 

This may come from the full analytic form of the remainder function at two loops i?g . The 
function Rq was calculated by Drummond, Hcnn, Korchemsky and Sokatchev [29] using the 
duality between the light-like Wilson loops and the MHV amplitudes, and then greatly simpli- 
fied by Goncharov, Spradlin, Vergu and Volovich [22] using the integral representation of Del 
Duca, Duhr and Smirnov [20j [21] . In the next section we perform the analytic continuation of 
the two-loop remainder function calculated by Goncharov, Spradlin, Vergu and Volovich to the 
region where ui — |ui |e — l27r , which corresponds to the Mandelstam channel in Fig. \Tjp in the 
multi-Regge kinematics. 

3 Analytic continuation 

In this section we discuss some details of the analytic continuation to the Mandelstam channel 
in the Regge kinematics. We also show how the kinematics determines the physical region of the 
cross ratios and establish the match between our picture and the one drawn by Alday, Gaiotto 
and Maldacena [30) . 

The result of Goncharov, Spradlin, Vergu and Volovich |22[ for the two-loop remainder 
function reads 



R, 



(2), 



(> 



3 ( 1 

(u!,u 2 ,u 3 ) = ^2 [ L ^( x i) x i ) - 2 Lizl( ' 1 ~ 



I (^2 Lia(l ~ V«<)) + Y A + {J 2 + C(2)) , (13) 



1 When the present manuscript was already at the last stage of the preparation, a new version of ref. [22) appeared. 
The non-analytic term \ was eliminated in the new version. This fact does not affect our result so that here we use 
the initial version of the GSVV formula. 
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where 



± ± ± u l + u 2 + U 3 - 1 ± \fK 
X; =mx , x = , (14) 



and A = (iti + u 2 + U3 — l) 2 — 4miM2U3. 
The function L^(x + , x~) is defined by 



3 (—\) m 1 
L 4 (z+z-)= £ ^J^log( 2 ;+x-r(4- m (x + )+4- m (^)) + ^ylog(x + ^) 4 , (15) 

together with 

£ n (x) = X - (Li„(x) - (-l)"Li„(l/x)) , (16) 

as well as the quantities 

j = ^(4(4)-4(*r)), (17) 

i=l 

and 

f -2 A > and u x + n 2 + U3 > 1, 
X = % 

+ 1 otherwise. 



(18) 



The result of the analytic continuation to the Mandelstam channel illustrated in Fig. [TJ), 
where u\ — |ui|e~' i27r , was presented by the us in ref. [25] and reads (see Appendices iDl and El 
for more details) 



+ y ln(|z| 2 |l - z\ 2 ) (lnzln(l - z) + lnz*ln(l - z*) - 2( 2 ) 

?7T ll — z\ 2 

+^ In l - T1 ^- (U 2 (z) + Li 2 (z*) - Li 2 (l - z) - Li 2 (l - z*)) 
2 \z\ 2 

+i2n (U 3 (z) + U 3 (z*) + Li 3 (l - z) + Li 3 (l - z*) - 2Q 3 ) . (19) 
In Eq. [19] we introduced complex variables 




to remove some square roots in the arguments of the polylogarithms (see Eq. I14j) . The phases 
4>2 and <p3 can be easily expressed in terms of the cross ratios m and have a meaning of the 
angles of the "unitarity" triangle illustrated in Fig. [2j More details on this parametrization are 
presented in the appendix [El 

The first term on RHS of Eq. [19] reproduces the leading logarithm term found by Bartels, 
Sabio Vera and one of the authors [7] in the BFKL approach as explained below. It is easy to 
see from Eq. [2D] that 

ln|z| 2 ln|l - z\ 2 = \nu 2 \nu 3 (21) 

and thus the first term on RHS in Eq. [TH equals to LLA in Eq. [TU] 

Other terms correspond to the next-to-leading logarithmic approximation (NLLA) and they 
present a new result, which is yet to be calculated using the BFKL technique. This analysis 
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(0,0) 



(1,0) 



Figure 2: The "unitarity" triangle. 



shows an agreement between the conjectured duality between the light-like Wilson loops and 
the MHV amplitudes, and the BFKL approach at the leading logarithmic level. 

The remainder function of Eq. [T5] in this channel in the multi-Regge limit is pure imaginary 
and symmetric under the substitution z f> 1 - z, which corresponds to the target-projectile 
symmetry in Fig. Q] Eq. 1191 vanishes for z — \ 1 or z — > 0, when the momentum of one of the 
produced particles fcj in Fig. Q] goes to zero, in an accordance to the expectation that in the 
collinear limit the six-point amplitude reduces to the five-point amplitude. 

Another useful form of the remainder function in Eq. [T5] can be written as 



n(2) LLA+NLLA I , , _; 2 7r 1 l w | 2 \ /-> mm, |2 i 



1 



1 2 



(22) 



~ In M 2 In 2 |1 + w\ 2 - ^ In 3 |1 + w\ 2 + wrln \w\ 2 (U 2 (-w) + U 2 (-w*)) 
-i27T {Li 3 (- W ) + Li 3 (-w*)), 
where the complex variable w is expressed in terms of the reduced cross ratios of Eq. [TTJ as 

1-z B+ „ l-z* B- 



U 2 ' Z* U 2 



(23) 



for S ± defined in Eq. ID.4l bv 



B ± _ l-u 2 -uz± -u 2 - u 3 ) 2 - 4u 2 u 3 ^ 

The complete discussion on the details of the analytic continuation and the z and w rep- 
resentations of the remainder function is presented in the appendices A-E, and here we only 
want to emphasize some important points. The analytic continuation was performed under an 
assumption that u\ + u 2 + u 3 < 1 to avoid a difficulty related to non-analyticity of y in Eq. 1181 
We made sure that after the continuation the remainder function does not have any singularities 
on the border of this region and therefore it is valid in the whole physical region. 

The variables u 2 and u 3 in Eq. [TT] are also cross ratios in the transverse momentum space 
as can be seen from Eq. [12] defining the dual coordinates in the transverse momenta space as 
illustrated in Fig. [3] 

In terms of the dual coordinates the reduced cross ratios in Eq. [12] read 

\x 0B \ 2 \x 'a\ 2 - \x a\ 2 \x 'b\ 2 

" 2 = i 121 12' U3 = I 121 12' " b > 

\xab\ \xqqi\ \xab\ \Xoo>\ 

Due to the Mobius invariance we can put 

xa = 1, xb = 0, xqi — oo, xq = z, (26) 
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Figure 3: The dual coordinates of the transverse momenta. 



then 



U2 



t*3 = |l 



(27) 



for z given by Eq. [20] This imposes a restriction on the possible values of the reduced cross 
ratios as illustrated in Fig. |4] (see ref. [25] for more details) . 



~ 2 

U 2 



N 

e x 



~ 2 

U 3 



Figure 4: The region of the reduced cross ratios where the analytic continuation is performed. 

The region A in Fig. [4] is the region of possible values of the reduced cross ratios, which 
correspond to the particle momentum parametrization. Its subregion B is the region, where the 
analytic continuation is performed. As it was already mentioned we made sure that our result 
is valid in the whole region A. 

The same conclusion concerning the physical region of the cross ratios can be reached 
using the parametrization of the cross ratios introduced by Goncharov, Spradlin, Vergu and 
Volovich |22) . One can parametrize the cross ratios by six complex variables, namely 

^23^56 ^16^34 212^45 , na \ 

Ul = , U 2 = , U 3 = , (28) 

^25^36 ^14^36 Z14Z25 

where = Zj — Zj . In this parametrization the square roots in the arguments of the remainder 
function Eq. Q2] disappear and xf are rational functions. Exploiting the conformal invariance 
we can set 

24 = 0, 25 = 1, zq — 00. (29) 
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Then Eq. [28] reads 



z 3 - z 2 z 3 z\ - z 2 
ui = , u 2 = — , u 3 = — -. 

l-2 2 Z 2 Zi(l-Z 2 ) 



(30) 



Solving Eq. [30] for Zi we obtain some square roots that determine the physical region. For 
example, one of the solutions is given by 



-1 + ui + u 2 + u 3 ± - ui - u 2 - u 3 ) 2 - 4«i 



u 2 u 3 



2u 2 u 3 



(31) 



The argument of the square root coincides with A = (1 — u\ — u 2 — u 3 ) 2 — Auiu 2 u 3 defined in 
Eq. [TJ] The surface A = determines the boundary of the space of the physical values of the 
cross ratios. This surface is depicted in Fig. [5] 




Figure 5: The surface A = is the boundary of the physical values of cross ratios. 



The same plot, but for AdS^ surface was obtained by Alday, Gaiotto and Maldacena [50] 
introducing momenta parametrization of the cross ratios. The space of the physical cross ratios 
in the unit cube is inside the "bag" in Fig. [5] To find a match between our picture of the physical 
region depicted in Fig. [4] and the "bag" in Fig. [5] we draw the surface A = in the coordinates 
1 — iti, Vva and ^/ui} as illustrated in Fig. [6] 

It is clear from Fig. [5^ that in the Regge limit, when (ki+k2) 2 /s2 — 1—ui — > the boundary 
of the surface becomes a semi-infinite strip in accordance with Fig. [4] If one relaxes the Regge 
kinematics this region becomes a closed cigar shaped region as follows from the geodesies in 
Fig. Ek. 

In the next section we calculate the next-to-leading-order impact factor from our result of 
the analytic continuation of the remainder function in Eq. 1191 
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(a) (b) 

Figure 6: The surface A = in the coordinates 1 — u±, \fu^ and \fu^. In the Regge limit, when 
(ki + k2) 2 /s2 — 1 — u\ — > we obtain the region A in Fig. [U Once the Regge kinematics is relaxed 
we get a physical region of a cigar shape instead of the semi-infinite strip in Fig. [H The figures (a) 
and (b) depicts a different view of the same surface A = 0. 

4 NLO impact factor 

In this section we calculate the next-to-leading (NLO) impact factor appearing in the BFKL 
approach. We begin with Eq. [4] and Eq. [5j which define the imaginary part of the BDS violating 
term with logarithmic accuracy. We are interested in generalizing Eq. [5] to include also next- 
to-leading in logarithm of the energy (NLLA) terms. Taking into account NLLA corrections 
corresponds to relaxing multi-Regge kinematics to quasi-multi-Regge kinematics (QMRK) for 
the intermediate particles in the unitarity relation for the amplitude. 

The next-to-leading corrections to Eq. [5] are of two kinds: the NLO corrections to the 
impact factors of the BFKL ladder and to the kernel of the BFKL equation. The later was 
calculated in QCD by Fadin and Fiore [311 I3"2] and in M = 4 SYM it can be extracted from 
QCD calculations applying a principle of the maximal transcendentality proposed by Kotikov 
and one of the authors |33j . The maximal transcendentality principle was successfully used, 
for example, to predict the anomalous dimension up to six loops in J\f = 4 SYM [33J [351 123 
l37l l38l l39l l40l . Another sort of NLO corrections to Eq. [5jis the corrections to the Reggeon- 
Reggeon-Particles (RRP) impact factor of the BFKL ladder, which were never calculated before 
for the octet channel It worth mentioning that there is some ambiguity in the higher order 
terms, namely, the corrections can be redistributed between the impact factor and the BFKL 
Kernel, which does not affect the form of the amplitude provided the corrections are included 
in a consistent way. 

We are interested in the NLO contribution to the impact factor. We write Eq. [5J for the 
leading logarithmic contribution as 

2 The NLO corrections to the impact factor in the singlet channel were found by Balitsky and Chirilli |41ll42j . 
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A 



2-S-4 



= 1 E (-D 



dv 



( qzh 



k 2llJ \k2qi, 

°° />oo 

I 2 / dv (2 Xl LLA K2 X 2 LLA ) ((1 - u,)-^ - l) , (32) 



n— — oc 



where the leading-log impact factors are given by 



1 



1 



2 (ii/. 



X2 o 



2 (--f) 



Oi 

as 



Si 
kt 

iv-\- 



(33) 



The functions Xi LA an( i X% LA are a convolution of the octet BFKL eigenfunction and the 
corresponding impact factor, but for the purpose of our discussion we call them impact factors in 
the v, n representation. For more details regarding the rigorous definition of the impact factors 
the reader is referred to ref. [7] . The factor of two accompanying X f LA m Eq. [32] is introduced 
to match the notation in Eqs. 90-93 of ref. [7]. 

In the appendix F we found that the NLLA term of the remainder function at two loops in 
Eq. [22] can be written as 



Rf ] NLLA ( Me"^, -f-^, ) (34) 




This allows to modify the leading logarithmic expression Eq. [32] to account for the next-to- 
leading in In S2 (NLLA) corrections at all orders of the perturbation theory. 

Before we proceed there is one fine point to be clarified. According to the factorization 
hypothesis |14) the all-order remainder function R is defined by 

M 2 _> 4 - Mg^fR (35) 

It was argued by one of the authors [28] that provided the factorization hypothesis holds, a 
significant information about the remainder function can be obtained from the analytic prop- 
erties of the BDS formula. In particular, in the region under consideration, where s, S2 > 
and S\, S3, S012, S123 < in the Regge kinematics the remainder function at all orders of the 
perturbation theory is given by the dispersion relation (se^l Eq. 50 of ref. [28] ) 

/ioo j 
-/( W )e-~(1- Mi r, (36) 
-ioo 27r * 

where the phases 5 and uj a b represent contribution of Regge poles obtained directly from the 
BDS formula and are given by 

= — in(u 2 u 3 ) , u ab = — In — . (37) 
8 8 1*2 



3 The function ii is denoted by c, and the reduced cross ration ui and U3 are $2 and 4>3 i n the notation of ref. |28] 
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The coefficient is the cusp anomalous dimension known to any order of the perturbation 
theory. The only unknown piece in Eq. [36] is the real function f(u), which has the Mandelstam 
cut in u and depends only on the transverse momenta and has no energy dependence. In the 
leading logarithmic approximation f(w) is given by 



f ^)= a - v u \ , t^my-u^L) - ( 38) 



9 



... dv (2 X i LA )(2 X ^ LA ), 



where oj(v, n) is defined in Eq. [5] Therefore in RHS of Eq. [35] the integral over uj gives correctly 
the leading asymptotics of imaginary and real parts of the amplitudes (see Eq. [SJ. 

The expression for J lla (uj) is read out from Eq. [33] and can be generalized to include 
subleading contributions. By analogy with Eq. [38] we write 

n— — oo v y 

where both the impact factors Xi an d the BFKL energy u)(v, n) include the next-to-leading 
corrections. They are defined by 

Xi=X? LA + X? L ° (40) 

and 

Q(u,n) = u)(v,n)+u NLO {v,n). (41) 

The expression for lo nlo {v 1 n) can be found from the next-to- leading corrections to the octet 
BFKL Kernel calculated by Fadin and Fiore [31] [32] , and the missing NLO impact factors we 
can readily extract from Eq. 1341 In the appendix F we find that the next-to-leading impact 
factors Xi are given by 



nlo _ " 1 ( 1 n \ I 9i \ 2 ( 5i 




— IV — 4* / *\ — iu+-$ 



V. =7" ri^„- T - -|( (-|) <42i 



and 

X2 ""i(^r-""S(^fJ UJ ' (43) 

where E v<n is defined in Eq. [5] An important feature of Xf LO is that both Eq. [4j and Eq. [33] 
do not have holomorphic separability, namely either w + n/2 or — ii/ + n/2 cannot be assigned 
only to one of the projectiles. It is worth emphasizing that the NLO impact factors x^ are 
factorized in the product of the Born impact factors in Eq. [33] and a term expressed through the 
eigenvalue E V}Tl of the BFKL equation in LLA. The form of the NLO impact factor in the v, n 
representation resembles the three-loop remainder function in LLA, emphasizing the intimate 
relation between the two as discussed in the next section. 



5 Three loops in LLA and NLLA 

In this section we calculate the three-loop leading logarithmic (LLA) contribution to the re- 
mainder function of the six-point MHV amplitude and find the real part of the subleading in 
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lns 2 /so (NLLA) term. In the leading logarithm approximation (LLA) one neglects all terms, 
which are not enhanced by the logarithm of energy. In our case the main contribution comes 
from the Mandelstam cut in the variable U) canonically conjugated to In s 2 (see Fig. [T]) and the 
leading terms are those which have each power of the coupling constant accompanied by the 
same power of the logarithm of the energy lns 2 /so or equivalently by ln(l — u\), because in the 
Regge limit 1 — U\ ~ (kx + k 2 ) 2 /s 2 . The all-loop LLA contribution to the remainder function 
in the Mandelstam channel (see Fig. [TJj) is given by Eq. @] and Eq. [5j The three-loop term is 
obtained expanding s 2 ^ v,n ^ in Eq. [5] in powers of the coupling a and reads 

{3)LLA iA 2 % 4 i 2 ^ /■ (-!)" (m\ iV -* {9skx\ iu+ f A ,, 



n— — oo 

oo 



1 1 J/l „. A f J.. ( X ) L-2 



= - 1^(1 - «x) 2J / ^-jf^^n M 7 , 

n=— oo ^ 4 

where the complex variable u; is defined by 

W = , W = -t 3 - . 

k 2 qi k^ql 

The expression in Eq. [44] is calculated in the appendix [G] and has the form 



(45) 



^(3) LLA = lA (^J a 3 = i7f l m 2 (1 _ ^ n | w |2 ^ |j_ + ^2 _ 2 ^3 ^ + ^2 (4g) 

-i In 2 |u;| 2 ln |1 + w| 2 + - ln \w\ 2 (Li 2 (-w) + Li 2 (-tfl*)) - U 3 (-w) - Li 3 (-w* 

The LLA three- loop remainder function in Eq. 01] is pure imaginary and vanishes at w — > 
and w; — )■ oo. It is invariant under w — > 1/w transformation, which corresponds to the target- 
projectile symmetry. 

The next-to-leading in the logarithm of the energy ln s 2 /s — — ki(l — Ui) contribution can 
be obtained from Eq. [35] and Eq. 1391 Unfortunately we do not have an explicit expression for 
u)(v, n) beyond the leading order, which is necessary for this calculation. However it turns out 
that we can find asymptotic behavior of the real part of the NLLA remainder function at three 
loops, because it does not require any knowledge of the higher order corrections to u}(y,n). 
?St [R^ NLLA ^j is calculated expanding Eq. (35] in powers of a. The details of this analysis are 
presented in the appendix [HI and the result is given by 

K(4 3) NLLA ) = ^ ln(l - Ul ) (in \w\ 2 In 2 |1 + w\ 2 _ |ln 3 |1 + w\ 2 (47) 

-i In 2 H 2 ln |1 + w\ 2 - ln |w| 2 (Li 2 (-w) + Li 2 (-w*)) + 2Li 3 (-iu) + 2Li 3 (-u>*) 

Note that SR(i4 3 } NLLA ) resembles very much the form of R^ LLA in Eq. lG~18l as one could 
expect from Eq. IH.4I The complex variables w is expressed in terms of the reduced cross ratios 
u 2 = u 2 /(l — Mi) and u 3 — 1(3/(1 — ui) as 

w = = — , w* = — = — (48) 

Z U 2 Z* U 2 

for S ± defined in Eq. ID.4l bv 



, 1 - U 2 - U 3 ± - U 2 - U 3 ) 2 - 4U 2 U 3 
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6 Conclusion 



In this paper we present some details of the analytic continuation of the GSVV [32] formula 

(2) 

for the two-loop remainder function i?g to a physical region of the 2 — > 4 particle scattering, 

(2) 

where i?g gives a non-vanishing contribution in the Reggc limit. Wc hnd that after the analytic 
continuation the remainder function reproduces the leading logarithmic (LLA) term calculated 
by Bartels, Sabio Vera and one of the authors [7] in the BFKL approach. We also find a 
term subleading in the logarithm of the energy (NLLA) and extract the next-to- leading (NLO) 
impact factor used in the BFKL technique. The BFKL approach allows to calculate the LLA 
contribution to the remainder function at any order of the perturbation theory. The three loop 
LLA remainder function as well as the real part of the three-loop NLLA remainder function are 
calculated in the BFKL technique and presented in this study. We find that the all-loop LLA 
and the two- loop NLLA terms of the remainder function are pure imaginary in the Regge limit, 
while starting at three loops the NLLA remainder function develops a non-vanishing real part 
in the Regge limit after the analytic continuation to the relevant physical region. 
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A Polylogarithms of U{ 

We perform an analytic continuation of the GSVV expression Eq. [Tj|]in the variable u\. The 
functions which do not depend on u\ remain the same after the continuation. In the multi Regge 
kinematics given by Eq. [5J they can be simplified as follows 

Li 2 (l- — J ~-C2-Jln 2 (l-ui)-ln(l- Ul )ln U2 -iln 2 u 2 (A.l) 
V U2 J 2 2 

and 

Li 4 (l - ^) - - | ln2 (! - «i) - YA 1r4(1 " Ul) ~ Ca ln(1 " Ul) ln " 2 (A ' 2) 

~ ln 3 (l - ui) \nu 2 - y m2 «2 — ^ ln 2 (l ~ «i) hi 2 m 2 - ^ ln(l - u x ) In 3 u 2 - ^ In 4 u 2 . 

in terms of the reduced cross ratios defined in Eq. 1111 

The expressions for Li 2 ^1 — and Li 4 ^1 — are obtained from Eq. IA.1I and Eq. IA.2I 
replacing m 2 by M3 in the argument. 

The polylogarithms Li„ ^1 — should be analytically continued and in the multi Regge 
kinematics are given by 

— ln(l - i) ~ -i27r / — = -i27rln 1 - - — r (A.3) 

t Ji t v Kl/ 

= -i27r(ln(l - mi) + in) = -i27rln(l - m) + 2tt 2 . 

Note that we assign a phase ln(— 1) = +iir since the argument 1 — goes counterclockwise 
around the origin as we continue from ui = e z0 |i*i | to u\ = e _l27r |Mi| through u\ = e~ l7T \u\\ = 
— \ux\. In a similar way we find 







K) 


-I 







4 

— + i7 r 3 ln(l - mi) + 7T 2 ln 2 (l - ui) - — ln 3 (l - in). 
3 3 

± — 1 ™± 



B Polylogarithms of x ± and x 

The variables x 1 * 1 are defined as follows 



± Ml + M 2 + M 3 - 1 ± y/A ± ± 

£ = , x\ = UiX , (B.l) 

2mim 2 m 3 

where A = (ui + m 2 + u% — l) 2 — 4mim 2 M3. 

First we consider the logarithm of x~ /x + . This variable goes from the second quadrant of 
the complex plane (negative real and positive imaginary part) in the counterclockwise direction 
during the continuation crossing the real axis for the large negative values of the argument, then 
it again crosses the real axis for small negative values of the variable. Thus the argument crosses 
the branch cut twice and lnx~ /x + remains on the same Riemann sheet and does not acquire 
any imaginary part after the continuation. Namely, 



, x \x I 1 - U2 - U3 + V(I - m 2 - m 3 ) 2 - 4m 2 m 3 

m— =ln— — ~ , (1>.2) 

x F I 1 - m 2 — u 3 — — u 2 - m 3 ) 2 - 4m 2 M3 
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where the notation \x^\ denotes the fact that u\ in x^ is replaced by |ui|. 

For our purposes it is useful to note that ln(l — z) has the same cut structure as Li n (z). This 
can be seen from the series representation 

k=l 

and thus 

OO 

u 1 ( Z ) = j2y = - 1 < 1 -^ (R4) 

k=l 

B.l Logarithms of xfx~ 

The variables xf are defined in Eq. IB. II The functions of these variable that are present in the 
result of Goncharov et al. are ~Li n (xf) and Li n (l/xf) (n — 1...4). 

The logarithm of the product x\ ' x~ can be easily analytically continued noting that 

x t x i = . * = 1-3. (B.5) 

U l+ iU l+2 

We readily calculate 

\nxfxi = —i2ir — ln{*2 — lnu3 — 2 ln(l — \ui\), (B-6) 

\nx2X2 — i2ir + ln{t2 — ln{t3, (B.7) 

I11X3X3 — i2ir — ln-5,2 + hiU3, (B.8) 
because lnui ~ — i2ir after the continuation in the limit Eq. [21 

B.2 Poly logarithms of xf 

We start with xf . During the analytic continuation xf goes from the second quadrant of the 
complex plane (negative real and positive imaginary part) in the clockwise direction and then 
crosses the real axis between and 1. Thus Li n (xf) are not changed after continuation and can 
be simplified as follows. 

First we simplify the argument in the limit Eq. [2] separating the "longitudinal" ( 1 — Ui) and 
the "transverse" (v,2 and u 3 ) cross ratios. We write 



Af 



where is a function of only 112 and u 3 and is given by 



(B.9) 



A ± _ 1 - u 2 - u 3 =f y/(l - v.2 ~ u 3 ) 2 - 4m 2 % ^ iq 
1 2-S2U3 

Using this notation we can expand the polylogarithms as follows 

Lh(xf ) = Lixd.T+1) = - ln(l - |4|) ~ - InAf + ln(l - m), (B.ll) 
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Li 2 (4) = Li 2 (|a;+|) - - y - ^ In 2 A+ + In A+ ln(l - u x ) - ~ ln 2 (l - u x ), (B.12) 

Li 3 (x+) = Li 3 (|4|) ~ ~ In A+ - 1 In 3 A+ + ^- ln(l - «i) + | In 2 A+ ln(l - u#.13) 
-ilnA+ln 2 (l- Ul ) + iln 3 (l- Ul ), 

Li 4 (4) = LuGs+l) ~ - ^ In 2 - 1 In 4 A+ + ^ In A+ ln(l - u x ) (B.14) 

+ 1 In 3 A+ ln(l - Ul ) - 1L ln 2 (l - Ul ) - ~ In 2 A+ ln 2 (l - Ul ) 

+ ilnA+ln 3 (l- Ul )-^ln 4 (l- Ul )- 

The variable x^ also does not cross the branch cut of the polylogarithm (from 1 to +oo) 
since during the continuation it goes from the third quadrant of the complex plane (both the 
real and the imaginary parts are negative) in the clockwise direction and then crosses the real 
axis for negative values (never crosses the imaginary axis). The expansion of the poly logarithms 
of Xi is obtained from that of x^ replacing by . 

The polylogarithms of 1 /xf are analytically continued since the variable goes from the third 
quadrant of the complex plane in the counterclockwise direction, crosses the imaginary axis and 
then the real axis behind 1 when we go from u\ = e l0 |ui| to u\ = er %2 ^\u\\. The argument 
crosses the branch cut of Li n (l/x x ), and the direction of the rotation determines the sign of the 
phase of ln(l — l/x x ) as positive (+i2n). 

Using this notation we can expand the polylogarithms as follows 

Lii (-^J ~ -«2tt, (B.15) 

u > te) - f f - " ? - «> - ^ * ™ 

~ -i27r(ln(l - U X ) - hxAf + in) = -i27rrn(l -u x ) +«27rlnA+ + 2tt 2 . 

In the second line of Eq. IB.16I we used the fact that ln(— 1) = because the argument l/x x 
rotates in the counterclockwise direction around the origin. In a similar way we continue the 
rest of the polylogarithms. In general we can write 

Li„ (J^j ~ -i2n ^ l _ 1); (ln(l - u x ) - In A+ + itt)^ 1 . (B.17) 

The polylogarithms of l/a;J~ are not continued because the argument goes in the counter- 
clockwise direction in the left complex semi plane and never crosses the imaginary axis and thus 
also the branch cut of the polylogarithms. In our limit l/x x —¥ and thus 

Li„ (ij ~ 0. (B.18) 

From Eq. IB.15I and Eq. IB. 181 we calculate 

Li i ( -V) + Li i (— I ~ ~ i27T > ( B ' 19 ) 
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+L "(^) 







which can be checked by the direct calculation eliminating the square roots before the analytic 
continuation 

. u ( ! _ 1 ) _ i n ( X _ ±) = _ m (l-"2)(l-« 3 ) „ _ i2?r (B 20) 

This shows that we fix correctly the phase of the argument according to the direction of its 
rotation during the analytic continuation despite the presence of the square roots. 

B.3 Polylogarithms of xf 

The variable x 2 rotates from the second quadrant of the complex plane (negative real and 
positive imaginary parts) in the clockwise direction during the continuation, crossing the real axis 
in the left semi plane and never crossing the imaginary axis. This means that the polylogarithms 
of X2 remain the same after the analytic continuation since its argument never crosses the branch 
cut of the polylogarithm (from real values of 1 to +00). 

On contrary, x% rotates from the third quadrant in the counterclockwise direction crosses 
the imaginary axis and then the real axis beyond the value of 1, crossing the branch cut of the 
polylogarithms. Thus the polylogarithms of this argument are to be analytically continued. 

We start with simplifying Li„(a;^). For the sake of convenience, similarly to the previous 
discussion we write it as 

\xf\=Af=-Afu 2 , (B.21) 

where 



A ± _ U 2 +U 3 -l± y/(l -U 2 - U 3 ) 2 ~ 4M 2 % 

2 2u 3 ' 1 ' ' 

and At was defined in Eq. lB~Tol We use redundant definitions , which are expressible 
through each other solely for the sake of transparency of the intermediate calculations. The 
final result will be simplified using the relations between them. 

As we have already mentioned the polylogarithms of x^ remain the same after the continu- 
ation and are given by 

Li„ (x+) = Li n ~ Li„ (A+) . (B.23) 

The polylogarithms of x^ are to be continued as follows 

Li! (xj ) = - ln(l - £3 ) = - ln (! - A 2 ) - i2n (B.24) 

and 

r x 2 dt 

U 2 (x 2 )=- — ln(l - t) ~ Li 2 (A 2 ) - i2nkiA 2 . (B.25) 
Jo * 

In general 

Li„ (x 2 ) ~ Li„ (At) - z2tt * In"" 1 A£. (B.26) 

(n 1). 

Using Eq. lB~23l and Eq. lB~24l we write 

Lii (x£) + Lii (x 2 ) ~ - ln(l - A$) - ln(l - ) - i2?r = — z2tt + lnw 3 , (B.27) 
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where we used the identity 



(l-A+)(l-A 2 -) = i. (B.28) 

On the other hand we can simplify Lii (x^ ) +Lii (x 2 ) before the continuation and continue 
it after that 

Lii (x+) +Lii (xo) = -ln(l -xt) - ]n(l-Xo) = _ m ( 1 ~ ~ "s) _ _ i2 vr + lnu 3 . (B.29) 

U1U3 

This confirms our choice of the phase sign according to the direction of the rotation of the 
argument. 

The variable l/x 2 goes from the third quadrant of the complex plane in the clockwise 
direction during the continuation. It crosses the real axis in the left complex semi plane and 
thus never crosses the branch cut of the polylogarithms. The argument l/x 2 goes from the 
second quadrant of the complex plane in the clockwise direction during the continuation. It 
crosses the imaginary axis and the real axis, but never reaches the branch cut in our limit. Thus 
all polylogarithms of l/x 2 remain the same after the continuation and can be written as follows 

Li„ (X) ~ Li„ (-L) ■ (B.30) 



A 2 



As an example we calculate 



L « ) + Li ' {-k) - ■» (' - it) - " (' - -k) = < B - 31 > 

We can check directly the validity of this result by eliminating the square roots before the 
continuation 

B.4 Polylogarithms of xf 

The polylogarithms of x^ are analytically continued exactly in the same way as corresponding 
polylogarithms of x^. This can be shown by numerical calculations as well as explained on 
general grounds by the symmetry u 2 U3, which corresponds to target-projectile symmetry of 
the scattering amplitude. 

The polylogarithms of xf are obtained from that of xf by making a change U2 & "3 as well 
as replacing Af by Af, where 

Af = ^-Af = -Afu 3 . (B.33) 

B.5 Continuation of J 

The function J is defined in Eq. [T7] through the sum of the polylogarithm of xf . Using this 
definition we can readily find 

H M ^ + * V 4 ln If + * V " (B.34) 
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As it was anticipated in ref. [22 this contradicts the result of the continuation of the function 

ta£-h£} (B.35, 



found in the previous section (see Eq. IB.2[) due to the difference in the cut structure between the 
two. The correct analytic continuation was done using the definition of J in terms of Lii(x i ) 
in Eq. [TTJ 



C Leading Logarithmic Approximation (LLA) 

We want to extract the leading logarithmic term from the expression of the remainder function 
after the continuation. The leading logarithm of the energy In S2 is related to ln(l — u\) through 

|ki+k 2 | 2 

1 — «i oc (CI) 

The reduced cross ratios U2 and U3 depend only on the transverse components of the momenta 
of the external particles and are of the order of unity. It should be emphasized that only terms 
of the order of ln(l — u\) can contribute to the imaginary part of the remainder function at two 
loops. Higher order terms ln n (l — u\) (n > 1) would contradict the unitarity of the scattering 
matrix. These terms do appear at the intermediate steps of the calculations, but they all must 
cancel out in the final expression. 

From our previous discussions we see that only a few terms can have contributions to the 
LLA result 

\ 2 

1 

12 ' ' 




(C.2) 



We consider them in separate. For our purposes it is convenient to single out only LLA 
contributions produced in the process of the continuation. We use the fact that the remainder 
function vanishes before the continuation 

R ( q ] -> (C.3) 

in the limit u\ — > 1, 112 — > and M3 — > for U2 — U2/O- ~ ui) an d U3 = ^3/(1 — u\) being 
kept fixed and of the order of unity. We obtain the imaginary part of the remainder function 
in LLA by subtracting the contribution before the continuation from those obtained continuing 
the functions in the physical region where «i — > e~ t27r \ui\, keeping only terms accompanied by 
the power of ln(l — ui). We start with the first term in Eq. IC.2I 

where we used Eq. IA.4I and the fact that 

lim Li„ ( 1 - ) ~ 0. (C.5) 
In an analogous way we calculate the second term of Eq. IC.2I 

Mx+,xr) - M\x+\, K\)~ A ln(l - «i) - y ln 2 (l - «i) + f ln 3 (l - «x) (C.6) 

2 - 9 

— — ln(l — Mi) ln-U2 + — In (1 — ui) lnu2 + — ln(l — Ui) In U2 — ln(l — ui) I11U3 



In (1 — lii) lnu3 + — ln(l — Mi) ln{t2 lnit3 + — ln(l — lii) In M3. 
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The notation \xf\ denotes that fact that u\ is replaced by \u\\ in the argument. 
Finally, the last term in Eq. IC.2I is given by 




■In (1 — Mi) + — ln(l — u{) In u 2 — — In (1 — ui)\nu 2 — — ln(l — «i)hi u% 

+ — ln(l - Mi)lnu 3 - y ln2 ( 1 ~ u i) lnu 3 - -r ln(l - Ui)ln 2 u 3 . 

Adding up Eq. El Eq. \UM and Eq. [CJ] we get 

+ — ln(l - Mi)ln{i2lnit3. (C.8) 

This expression coincides with LLA term obtained by one of the authors j7j in the BFKL 
approach. 

D Next-to-leading logarithmic (NLLA) terms 



We have extracted the leading order term in the logarithm of the energy In s 2 — — ln(l— u\) of the 
imaginary part of the remainder function R e . The term we obtained from Eq. [IS] after analytic 
continuation to the physical region of u\ — > e~ ,27r |ui| reproduces the term calculated by one 
of the authors in the BFKL formalism. Unfortunately, due complexity of the calculations, the 
sub-leading terms in ln(l — u\) were not yet calculated in the BFKL approach. The comparison 
between the two approaches, the BFKL formalism and the Wilson Loop/Scattering Amplitude 
duality, is not full without matching the NLLA terms. In this section we calculated the NLLA 
terms from the analytically continued expression of Goncharov et al. given in Eq. Q2] that can 
be further confronted with the BFKL result once it is available. We follow the logic of the LLA 
calculations outlined above, and extract only the NLLA terms that appeared in the course of 
the continuation (subtracting the relevant values before they were analytically continued). As 
it was already mentioned this is possible to do because the remainder function vanishes in the 
limit of Eq. [5] before the analytic continuation. 

The expression of the remainder function is given in Eq. 1131 We calculate all contributions 
in separate leaving only the NLLA terms, i.e. those that are not accompanied by any power of 
ln(l-wi) 

lA / 1\ / 1 \nlo n 4 . , 

% — 1 v 7 L — l v 11/ 




L 4 (x+,xr) - u(\xt\, ldf=^ y + T ln2 B+ + J ln ' B+ - T" ln( " 2 " 3) (D ' 3) 

2 2 

-^lnB + ln(u 2 M 3 ) - ?ln 2 B + \a(u 2 u 3 ) -^-ln 2 (u 2 u 3 ) + ^ ln B + ln 2 (u 2 u 3 ) + ^ ln 3 (G 2 u 3 ), 
4 4 lb o 48 
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where we introduced 



B ± _ 1 - u 2 - u 3 ± ~ u 2 - m 3 ) 2 - 4m 2 m 3 p ^ 

and used its property 

= u 2 us (D.5) 

to eliminate B~ . 

Before we calculate the contributions from L^x^jX^) an( i L4 (a;^ , £3" ) we nno - ^ ne function 

J. 

1 B+ 

J~-ln—+m. (D.6) 

The expression in Eq. ID. 61 depends only on the " transverse" cross ratios u 2 and u 3 as one can 
see from Eq. ID. 41 The function \ defined in Eq. [18] has the same value x = 1 in all points on 
the circle u\ = \u\\e 1 ^ for |ui| + u 2 + 1*3 < 1 and thus does not posses any additional terms in 
the analytic continuation when the phase <f> changes from to —z27r. Now we can readily find 
the contribution from all terms that include J 

J l + x ±j* _ ( J l + X ±,A 1 „ _ £_ ^ B+ + a ln 3 B+ (D . 7) 

24 A 12 V 24 12 / 1 1-1 1 24 4 6 V ; 

2 ■ 2 

+ — lnB + ln(M 2 u 3 ) - ^j-ln 2 ^ + ln(w 2 u 3 ) - y^ln 2 (w 2 u 3 ) 



Note that 



-— lnB+ ln 2 (u 2 u 3 ) - — ln d (M 2 u 3 ). 



- ( X777C2 ) |« 1 =|ui| = 0. ! D. 



12" V 12 

Summing up Eq. ID.lllD.3| Eq. ID.7I and Eq. ID.8I we obtain a compact expression 

- ^ + I In 3 + £ In 2 f ^ - ln(«) - | In 2 5+ ln(«) + ^ lnB+ ln 2 (u 2 « 3 )(D.9) 
12 3 8 \ u 3/ 6 2 4 

As a last step in our analysis we calculate the contributions of \ja(x 2 , x%) and L,4,(x~3 ,x$). 
Namely, 
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* 3 ■ 3 

L 4 (x+,x 2 -) - U{\xt\, \x 2 \)"~ J '— - y In 3 B+ - ln 3 (B+ + S 2 ) + ^- lnS 3 (D.10) 
+ y ln(S+ + it 2 ) lnw 3 - — In 2 "3 + y ln(S+ + u 2 ) In 2 u 3 - — In 3 u 3 + — ln(u 2 u 3 ) 

9 9 

+ — In B + ln(S 2 S 3 ) — ln(B + + S 2 ) ln(S 2 S 3 ) — lnS 3 ln(S 2 S 3 ) 

— — ki(B + +M 2 )lnM 3 ln({t 2 M 3 ) - In 2 S 3 ln(S 2 S 3 ) + — — ln 2 (S 2 S 3 ) - — In 5+ ln 2 (S 2 S 3 ) 

Z 4 lb 4 

VK -. , i „ , . 2 / ~ - \ lS7T _ 2 . _ _ . Z7T 3 / ~ ~ \ Z7T 3 . i _ , 

+ — ln(S T +?i 2 )ln (it 2 u 3 ) + -3- hiu 3 ln (u 2 u 3 )-— In (w 2 u 3 ) — \n[B + u 3 ) 

o o 10 

2 ■ 2 

+ y lnS 3 ln(S+ + u 3 ) + y In 2 S 3 ln(B + + S3) - y ln(S 2 S 3 ) ln(B+ + S 3 ) 



— In u 3 ln(u 2 u 3 ) ln(B+ + u 3 ) + — In (m 2 u 3 ) ln(B + + S 3 ) + — Li 2 - 
z o Z \ u 2 

. 1 ~ T . { B+\ iTi" n _ / B+\ n\. ( B+\ . , _ _. / B+ 

+i7rlnu 3 Li 2 - — ln(u 2 u 3 )Li 2 — — — Li 2 — -i7rlnu 3 Li2 

V u 2 J 2 \ u 2 / 2 \ u 3 / \ u 3 

for, _ , T . / fl+N . T . / S+\ . T . / B+ 

+— ln(M 2 U3)Li2 — — - «ttLi 3 — - Z7rLl 3 — — 

2 \ «3/ V u 2 J \ U 3 

The contribution of L4 (2^" , Xg" ) is readily obtained from Eq. ID. 101 bv changing variables 
S 2 -H> S 3 . 

Summing up Eq. ID.9I and Eq. ID. 101 (together with S2 <~> S3) we get the NLLA part of the 
remainder function after the analytic continuation 



-— lnB + - — In 3 B + - — In 3 S 2 - Mrln£+ lnS 2 In S3 - — In 3 S 3 + — ln(S 2 S 3 )(D.ll) 
3 3 3 3 6 

-5ln 2 J B+ln(S 2 S 3 ) + % ln 3 (S 2 S 3 ) - l n\n ( ^) Li 2 (~) +inln ( ^) Li 2 



3 





' B+\ 


(-C)-^ Li3 ( 






, u 3 J 



2 6 ' V"3/ \ S 2 / VS3/ \ S 3 
— i27rLi3 

The expression of Eq. ID.lll is pure imaginary in the limit S2,3 > and S2 + S3 < 1 despite 
the fact that it contains a square root in its arguments through B + defined in Eq. ID. 41 It is 
also symmetrical with respect to the exchange of S2 and S 3 . 

Adding to Eq. ID. Ill the Leading Order result calculated in the previous section and given 
by Eq. IC.8I we obtain the final result 

+ yln(l-u 1 )lnS 2 lnS 3 (D.12) 

in 3 , 1 27r o 1 Z7r 3 , , , , in , 3 in 3 , , 

— — In B+ - — In 3 B+ - — In 3 u 2 - in In B+ In u 2 In u 3 - — In 3 u 3 + — ln(u 2 w 3 

3 3 3 3 

in i 2 n+ , ,~ - ^«r. 3 . , /S 2 \ T . ( B+\ , . , fu 2 \ T . ( B 



In i3 + ln(S 2 S 3 ) + — In (S 2 S 3 ) -27rln — Li 2 — — +Z7rln — Li^ 
2 6 \u 3 J \ u 2 J \U3/ \ u 3 

( B+\ , ( B A 

-«27rLi3 — - — — z27rLi 3 ( — - 



u 2 J \ u 3 

An important remark is in order. The last expression was calculated in the region S2 + u 3 < 1 
and, in principle, should be analytically continued to any other physical region. The first term 
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+^ln(l — lnu2 I11U3, which corresponds to the Leading Logarithmic Approximation is a 
smooth function also outside the region u 2 + u 3 < 1 since it does not have any singularities on 
the boundary of the region. This is not obvious for the rest of the NLLA terms, where individual 
terms do have branch points on the boundary of u 2 + u 3 < 1. However the singularities are 
canceled in the sum as can be shown introducing back B~ 

R(\ui\e~ z2 ' K , u 2 (l - ui), u 3 (l - ui)) ~ +— ln(l - m) \nu 2 lnu 3 + — In 3 u 2 



— In 2 uo In 

2 \U 3 

f f B+ 
-?2tt Li 3 

V V «2 



«2 



Z7rln 
■ Lis ( — 



Hi 

B~ 
11 2 



Li 2 



B+ 

U 2 



Li 2 



B~ 
u 2 



(D.13) 



It can be easily seen from the series representation of the polylogarithms that all square roots 
in the argument cancel out, and the expression in Eq. ID. 131 is also valid for u 2 + u 3 > 1, but 
only in the region A of the multi Regge kinematics shown in Fig. |4j Eq. ID.131 is the main result 
of this study. 

E Rq in complex variables 

In this section we eliminate the square roots in the arguments of the remainder function of 
Eq. ID.131 introducing complex variables 



z = V«2e 102 , 1 - z = Vwse" 403 . 
It is useful to calculate cos (02 — 03 ) and sin (02 — 03 ) from 

\z + (1 - z)\ 2 = 1 = U 2 + U 3 + 2^/U 2 '^/U 3 'cOs((j) 2 - 03)- 

We readily find 

I-U2-U3 



COS(02 







(E.l) 



(E.2) 



(E.3) 



and 



as well as 



sin(0 2 - 03 ) 



isin(0 2 - 03 ) 



yj4u 2 u 3 - (I-U2- U3) 2 



a/(1 - u 2 - u 3 ) 2 - Ait 



2M3 



(E.4) 



(E.5) 



2y/u 2 y/u 3 

With the help of Eq. Elan Eq. EH the function B ± , defined in Eq. EH can be written as 

B ± = e T#2 e T#3 | z | |i _ z \ (E.6) 

and thus 

B+ 1 - z B~ 1 - z* B+ z* B~ z 



H-2 



U 2 Z* 



u 3 1 — z* ' u 3 1 — Z 



(E.7) 
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Using Eq. IE.7I and the identities between Life of different arguments we write the expression in 
Eq.EIlas follows 

R(\ Ul \e- i2 *, |z| 2 , |1 - z| 2 ) ~ y ln(l - ui) In \z\ 2 In |1 - z\ 2 (E.8) 

-i47rC 3 + y In |z| 2 11 - z\ 2 (In z ln(l - z) + In z* ln(l - z*) - 2C 2 ) 
i7r II - z\ 2 

+ V 111 ^TiT" W*) + Li 2( Z *) - Li 2(! - *) - Li 2(l - **)) 

2 |z|2 

+i27r (Li 3 (z) + Li 3 (z*) + Li 3 (l -z) + Li 3 (l - z*)) . 

From Eq. IE.8l we see that the square roots present in B disappear and the remainder function is 
manifestly pure imaginary in LLA. The target-projectile symmetry U2 O u 3 , which is z O 1 — z 
symmetry in terms of the variables Eq. IE. II is also obvious in Eq. IE.8I 

Because of the holomorphic factorization of the impact factor in Eq. [5] it is more natural to 
express the final answer in complex variables 

u; = — , w* = (E.9) 

z z* 

Noting that the reduced crossed ratios ui and {t 3 are related to the transverse momenta (see 
Eq. [T^]) we can write 

qsh ^e"^ 3 1 -z 

= —==—— = =w (E.10) 

and thus Eq. [5] reads 

- 2 Ej- 1 ) ^ l— J l— J ( s2 - 1 

= ? E / ~T~~^ r w*~ 9 H w+f - 1 • can) 



n— — oo 



Eq. IE. Ill is explicitly symmetric inzf>l-z(!i)f>l/tu). We recast Eq. IE. 81 in terms of the 
variables w and ui* as follows 



R Ittile 



1 



|l + w| 2 ' |l + wp 



ln(l-ui)ln|l 



: ln 



(E.12) 



+ — In |™| 2 In 2 |1 + w\ 2 In 3 |1 + w\ 2 + iTrln \w\ 2 (Li 2 (-w) + Li 2 (-w*)) 

2 3 

-i2n (Li 3 (-w) +Li 3 (-«;*))• 



F NLO impact factor 

We wish to calculate inverse Mellin and Fourier transforms of the next-to-leading contribution 
to the remainder function. The form of the direct transforms in the complex variable w can be 
read out from the last line of Eq. IE. Ill and is given by 



f(w,w*) 



E 



dvp 



2ii/ i(j>n 



f(y,n), (F.l) 
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where w = pe^. The inverse transform thus reads 

w?L dp2 L 



-| POO p27T 

f(y,n) = — - dp 2 d(pp-^- 2 e-^ n f{pA), (F.2) 



which can be written as 



2 



1 ~, 



f(y,n) = j^Tt I d*w(w)-™-?-i(w*)- w +*-*f(w,w*). (F.3) 



The integration in Eq. IF. 31 is performed on the two dimensional plane in Cartesian coordinates 
Wi and u; 2 defined by w = Wi +iw 2 . We start with the logarithmic terms appearing in Eq. IE.12l 
The relevant logarithms ln fc |1 + w\ 2 can be obtained by differentiation of the power function 
(|1 + w\ 2 ) a and thus it is convenient to consider 

d 2 w \w\ 2b w* n \l + w\ 2a (FA) 

2f>U i „„|2a+2n| 



|26.„* n 1 1 i „,.|2a 

" (2tt)2 r(o + l) dz*J a 11 1+1 |z=z=J 



We introduce the master integral 



fl(o, 6; «) = - / d 2 x|x| 2a |z - x\ 2 \ (F.5) 



7T 



which corresponds to the one loop diagram with z being a momentum of the external particles. 
This integral is found by using the well-known formula of the momentum integration 

d d k _ m r(d/2 - Ai)r(d/2 - A 2 ) r(Ax + A 2 - d/2) 

n tv\ \tV\ \tvj \ \ \ /„2^, +Ao-rf/2 ' \ r - D ) 



(k 2 )^(( q -k) 2 )^ r(Ai)r(A 2 )r(d-Ai- a 2 ) ( q 2)A 1+ A 2 -d/ 2 

and it reads 

r(i + a)r(i + b) r(-i-a-b) 

9[a ' ' Zj ~ r(-a)r(-6)r(2 + a + b) |*|2(-i-a-6) ■ ^ 
Using the identity r(x)F(l — x) = 7r/sin(7rx) this can be written as 

gMz) = r 2 (i + q)r 2 (i + b) sin^sin^ (1+a+b , 

yv ' ' ' r 2 (2 + a + 6) 7rsin7r(a + 6)' 1 V ; 



From Eq. IF.4I and Eq. IF.5I we find a general expression for the inverse transform of the 
logarithms 



_ (-1)" d k d m r 2 (l + a)r(l + &)r(l + b-n) sinTrosinTrfr , 
y(k,m){y,n)- ^ dakdhm r(2 + 6 + a)r(2 + a + 6-n) ^sin^(a + 6) la=0 '" = ^ 1+ " /2 ( j 



Applying Eq. IF.9l for k — 1 and m — we get (for n^O) 



2vrln|l + H 2 ^— (F.10) 

V "r 4 



in full agreement with the Born term of the six-point amplitude. 
For fc = 2 and to = we get from Eq. IF. 91 

2 7 rln 2 |l + W | 2 ^2l^^^ + ^+^('l-*^+^)-2^(l)), (F.FI.) 
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where we used the identity ij)(z) = ip(l — z) — Trcotnz. 
Plugging k = 3 and m = in Eq. IF. 91 we get 



2 



2.1n 3 |l + H 2 ^-j r ^(- 7 r 2 + 3(^ [iv + M) + ^ (l - iv + M ) - 2^ (1] 

+3 ( V' ( + y) - ( l - w + t) + 2 ^ ) (F ' 12) 

and for k — 1 and m — 1 one obtains 

27rln|w;| 2 ln|l + w| 2 => (-1)™ 2W _ . (F.13) 

(*» + f ) 

Therefore from Eq. IF. 131 and Eq. IF. Ill we get 

2 



27iTn|l + io r In 



1 

1 + - 



= 27rln 2 |1 + w\ 2 - 27rln|w;| 2 ln|l + w\ 2 =4> (F.14) 



in an agreement with the analysis of ref . [7] . 
In a similar way we calculate 

|2 i„2 



2 / V 2 



2 ^v^n 
T 



2ir\n\w\ z \n z |l + w|^=> (F.15) 

H 

2 



-V' ^1 - if + y J + ^' ( w + y 
The rest of the terms in Eq. IE. 121 arc found by noting that 

27rln|w;| 2 (Li 2 (-io) +Li 2 (-io*)) - 4tt (Li 3 (-u>) +Li 3 (-w*)) =► * — '- — r . (F.16) 

(- 2 + x) 

This can be easily checked by calculating residue at v = ±z|n|/2 and summing over n. Both 
of the terms on LHS of Eq. IF. 161 have poles double poles at v = for n ^ (for n = 
the remainder function vanishes), which correspond to the infrared divergencies absent in the 
remainder function. Due to the special coefficients of these terms appearing in Eq. IE. 121 the 
infrared divergency is canceled in the final expression. This fact suggests that the relative 
coefficients of the individual terms can be fixed demanding the absence of the poles at v = for 
n y^O, together with w — > l/w symmetry. 

Gathering together the inverse transform of all terms in Eq. IE. 121 we finally obtain 

*»" < M^*,™ - US- ( «, - . OUT) 



|1 + U'| 2, |1 + »IV 2 „2 + !Ei ^ "'" 4(„2 + ! £) 2 



where E UJl is given by Eq. [8] This expression vanishes for n = 0, which corresponds to absence 
of the infrared divergencies in the remainder function; it is symmetric in n —¥ —n and v — \ —i>, 
which implied by the target-projectile symmetry w — > l/w. 
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From the definition of the impact factors in v, n representation given by Eq. [33] we read out 
the form of the next-to-leading-order (NLO) impact factor 



NLO _ a 1 / -2 1 n W ?3 \ / 9 3 



In the next section we use Eq. IF. 171 to calculate the three loop leading-log contribution to 
the six-point amplitude. 

G Three loop contribution in LLA 

The general expression for the leading logarithmic contribution to the imaginary part of the 
remainder function at any number of loops is given by 

— 1 tj-»"C^h (ST (HH (*"'-) <*•> 



= i e y (™t~* H w+t ( s 2 ( " n) - 1) , 



4 

where the u>{v,n) is related to the eigenvalue of the octet BFKL Hamiltonian E vn by 

u>(v,n) = -aE Vtn (G.2) 

with 



2^ + nd \ 2 J \ 2 



The complex variable w is defined through 



w 



= — and w = — — . (G.4) 



According to the discussion presented in appendix E in the Regge limit w can be written as 



with 



ir = \l — e"^° (G.5) 

U2 



cos(0 2 - 03) = 1 - U / ^J l3 (G.6) 

2y / U 2 U 3 



and 



sm(02 - 0s) = ~ 7== (G.7) 
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for 

u 2 = , u 3 = . (G.8) 

1 — «i 1 — Ul 

The LLA three-loop contribution to the remainder function is obtained from Eq. IG.2I 
by expanding in powers of the coupling constant a as follows 

3 d( 3 ) LLA 3 oo ^oo , 

" 6 -- — = = t ln2 52 E / ~rb? ^*r~ f H w+t 



4 



=-cx> J -°° ^ + X 

ln 2 (l -«!) £ ("!)" / "ITT^ H W+ " 



4 

where we used the fact that in the Regge limit 1 — u\ ~ (ki + k2) 2 /s2. 

The integral in RHS of Eq. IG.9I can be easily obtained using the calculations of the previous 
section. Namely, we have shown that Eq. IF. 171 gives the v, n representation of the NLLA 
contribution to the remainder function at two loops, namely 



n— — oo 



TV 1^- 



4 



We calculate in separate the transform of the second term in the brackets in Eq. IG.10I 



n— — oo 



(- 2 + f)' 



The integral in Eq. IG.11I can be calculated using the residue theorem closing the contour either 
in the upper semiplane for poles v = i\n\/2 and multiplying the residue by i27r, or in the lower 
semiplane for poles v = — i|n|/2 and then multiplying the residue by — i2tt. The result has 
w <-> 1/w symmetry so that it is enough to consider only contributions for \w\ < 1. The residue 
at v = —i\n\/2 for \w\ < 1 reads 

3 (-l) n 7r(w*)l n l 3(~l) n 7r(w*)l"lln|w| 2 (-l) n 7r(>*)l n l In 2 |w| 2 



2 |n| 3 4n 2 8|n| 

The summation over n (for n > 0) is readily performed using the series representation of the 
polylogarithms Li n (x) = J2kLi x k /k n and we get 

-i2^Y Res [ ( " X) " - [ A (w*) lv -% (u>) itf+ * ,-i\n\/2 ) (G.13) 

h W + X 2 V 4 (- 2 + f)7 / 

1 3 3 

= — 7r In |iti| 2 ln(l + w*) H — 7r ln |w| 2 Li2(— w*) ttLi 3 (—w*). 

8 4 2 

The contribution from the sum over negative n is added by substitution w* —> w and we obtain 



n— — oo 

l 



^ 2 + x 2 V 4 (^ 2 , , 



7rln \w\ \n\l + w\ +-7rln|w| (Li 2 (— «;) - Li 2 (— «;*)) - -7rLi 3 (— to) - -nLi 3 (—w*). 
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From Eq. IG.9I it follows that for the three-loop LLA contribution we need to calculate the 
following expression 

£ / dv-t^El n (w*)*-* ( W f +§ ■ (G.15) 

n— — oo 4 

This can be obtained by subtracting Eq. IG. 141 from Eq. IG. 101 In the previous section we found 
that 

i t I {k. - kyyr^) or* {g.w) 

n— — oo ' 



+ f ) 

2^ (Li 3 (-tu)+Li 3 (-u;*)) 



= J In |w| 2 In 2 |1 + w| 2 - ^ In 3 |1 + w| 2 + tt In |w| 2 (Li 2 (-iu) + Li 2 (-w*)) 



and thus we write 



-J In 2 |u;| 2 In |1 + w\ 2 + J In | W | 2 (Li 2 (-«;) + Li 2 (-«;*)) - | (Lis(-«0 + Li 3 (-™*)) . 

Finally from Eq. IG.9I and Eq. IG.17I we obtain the remainder function at three loops in the 
leading logarithm approximation (LLA), namely 

a3R (3) LLA = jA (3^ 4 = ln 2 (1 _ A |w| 2 ^2 ^ + ^2 _ 2 ^3 ^ + ^2 (G Jg) 

-~ In 2 H 2 In |1 + w| 2 + ~ In \w\ 2 (U 2 (-w) + Li a (-u>*)) - Li3(-u;) - Li 3 (-w*) 
The complex variables w is expressed in terms of the reduced cross ratios of Eq. [TTJ as 

w= = — , w = — = — (G.19) 

z u 2 z* u 2 

for B ± defined in Eq. ID.4l by 



, l-u 2 -u 3 ± y/(l-U 2 - u 3 ) 2 - 4u 2 u 3 



B ± = ~< -o - V V- ^ ^ (Q 2Q) 



H The real part of the remainder function at three loops 
in NLLA 

In this section we calculate the real part of the remainder function at three loops in the next- 
to-leading logarithmic approximation (NLLA). The expression for 3t(Rg NLLA ) is obtained 
expanding the dispersion relation Eq. [35] in powers of the perturbation expansion parameter a. 
It is worth emphasizing that the calculation of $l(R^ NLLA ) does not require the knowledge of 
currently unavailable subleading corrections to the BFKL eigenvalue oj(v,n). We plug the LLA 
function f LLA (uj) of Eq. [38] in Eq. l36l and expand it in a to the third order 

farfa'jtf' + « 3 i?f - ^ * 2£(ln(l - + z-) 2 f] I dv-t^El^ (w*)^ , (H.l) 

n— — 00 ' 4 
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where w and E v ,n are given by Eq. IG.4I and Eq. [8] respectively. The phases 5 and ui a b of Eq. [37] 
can be written as 

o = - ln(n 2 ti 3 ) = - In -— — ^, w ab = - in — = - In \w\ H.2 
2 2 |1 + w| 4 2 u 2 2 

using the leading order term for the cusp anomalous dimension — 4a. The equation Eq. IH.ll 
is valid only for the LLA term and the real part of the NLLA term of the remainder function 
at three loops. Solving it for the LLA term we get 



p (3) lla _ i 



ln 2 (l- Ul ) E [dv\^(w*f-H™) iV+9 El,n (H.3) 



n— — oo 



in full agreement with Eq. IG.91 

Next we solve Eq. IH.ll for the real part of the NLLA remainder function 

. ID (2) UA oo „ / ,\ np 2 

K(i? ( 3) NLLA) = 1 H1 _ ui) ^2 / 4,i=^!t H *+* 

^ n— — oo ^ 4 

■ c D (2) LLA „ D (3) LLA 



i ln(l — i*i ) 



(H.4) 



From Eq. IH.4I we see that the next-to-leading logarithmic contribution is related to the lead- 
ing logarithmic terms at two and three loops. The function LLA was found using BFKL 
approach in ref. [7] and given by the first term on RHS of Eq. [22] 



R (2) i^ = !E ]n (i_ Ul )i I1 |i + 1i ,|a] n 



1 2 



(H.5) 



The LLA remainder function at three loops R^ LLA was calculated in the appendix [G] and is 
given by Eq. IG.18I Summing up all terms in Eq. IH.4I we readily obtain 



2 , 2 



$t(Rg> JV ^) = — ln(l - ui) ( In H ]n z \l + w\ z - - ln J |1 + w\ z (H.6) 



? (3)ArLL ^-^ln(l- Ml )(ln| W | 2 ln 2 ,. 

1 



-| In 2 \w\ 2 In |1 + w\ 2 - In |w| 2 (Li 2 (-w) + Li 2 (-w*)) + 2Li 3 (-w) + 2Li 3 (-w*)J . 

Note that 3?(i?^ 3) NLLA ) resembles very much the form of R^ LLA in Eq. IGT81 as one could 
expect from Eq. IH.4I The complex variables w is expressed in terms of the reduced cross ratios 
of Eq. [11] as 

l-z B+ ^ \-z* B- 

w* = — = — (H.7) 



Z U 2 Z" U 2 



for S ± defined in Eq. ID.4l bv 



D± 1-U2- U3 ± y/(l-U 2 - U3) 2 - 4u 2 U 3 / j j D ^ 

B = . (H.8) 



31 



References 

[1] S. J. Parke and T. R. Taylor, Phys. Rev. Lett. 56, 2459 (1986). 

[2] C. Anastasiou, Z. Bern, L. J. Dixon and D. A. Kosower, Phys. Rev. Lett. 91, 251602 (2003) 

[3] Z. Bern, L. J. Dixon and V. A. Smirnov, Phys. Rev. D 72, 085001 (2005) 

[4] L. F. Alday and J. Maldacena, JHEP 0711, 068 (2007) jarXiv: 071 0.10601 [hep-th]]. 

[5] J. Bartels, L. N. Lipatov and A. Sabio Vera, Phys. Rev. D 80, 045002 (2009) 

[6] O. Stcinmann, Helv. Physica Acta 33 (1960) 257, 349. 

[7] J. Bartels, L. N. Lipatov and A. Sabio Vera, Eur. Phys. J. C 65, 587 (2010) 

[8] L. N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338; 

V. S. Fadin, E. A. Kuraev, L. N. Lipatov, Phys. Lett. B 60 (1975) 50; 

E. A. Kuraev, L. N. Lipatov, V. S. Fadin, Sov. Phys. JETP 44 (1976) 443 ; 45 (1977) 199; 

I. I. Balitsky, L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

[9] J. Kwiecihski and M. Praszalowicz, Phys. Lett. B 94 (1980) 413. 

[10] J. Bartels, Z. Phys. C 60, 471 (1993). 

[11] L. N. Lipatov, J. Phys. A 42, 304020 (2009) 

[12] R. C. Brower, H. Nastase, H. J. Schnitzcr and C. I. Tan, Nucl. Phys. B 822, 301 (2009) 
larXiv:0809.l632l [hep-th]]. 

[13] J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, Nucl. Phys. B 826, 337 
(2010) 

[14] L. F. Alday and J. M. Maldacena, JHEP 0706, 064 (2007) |arXiv:0705.0303l [hep-th]]. 
[15] M. Drummond, G. P. Korchemsky and E. Sokatchev, Nucl. Phys. B 795, 385 (2008) 
[16] A. Brandhuber, P. Heslop and G. Travaglini, Nucl. Phys. B 794, 231 (2008) 
[17] N. Berkovits and J. Maldacena, JHEP 0809, 062 (2008) 

[18] J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, Phys. Lett. B 662, 456 
(2008) 

[19] J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, Nucl. Phys. B 815, 142 
(2009) 

[20] V. Del Duca, C. Duhr and V. A. Smirnov, JHEP 1003, 099 (2010) |arXiv:0911. 53321 [hep- 
ph]]. 

[21] V. Del Duca, C. Duhr and V. A. Smirnov, JHEP 1005, 084 (2010) jarXiv:1003. 17021 [hep- 
th]]. 

[22] A. B. Goncharov, M. Spradlin, C. Vergu and A. Volovich. la"rXiv:1006. 57031 [hep-th] . 

[23] V. Del Duca, C. Duhr and V. A. Smirnov, JHEP 1009, 015 (2010) |arXiv:1006.4127l [hep- 
th]]. 

[24] L. F. Aldav JarXiv:1009.m0l [hep-th]. 

[25] L. N. Lipatov and A. Prvgarin. larXiv:1008.10i"6l [hep-th]. 

[26] R. M. Schabinger, JHEP 0911, 108 (2009) 

[27] J. Bartels, J. Kotanski and V. Schomerus. larXiv:1009!3938l [hep-th] . 
[28] L. N. Lipatov. larXiv:1008. 10151 [hep-th]. 

[29] J. M. Drummond, J. Henn, G. P. Korchemsky and E. Sokatchev, Nucl. Phys. B 815, 142 
(2009) larXiv:0803. 14661 [hep-th]]. 



32 



[30] L. F. Alday, D. Gaiotto and J. Maldacena. larXiv:0911^4708l [hep-th]. 

[31] V. S. Fadin and R. Fio re, Phys. Lett. B 610, 61 (2005) [Erratum-ibid. B 621, 61 (2005)] 
|arXiv:hep-ph/0412386) . 

[32] V. S. Fadin and R. Fiore, Phys. Rev. D 72, 014018 (2005) j larXiv:hep-ph/0502045| . 

[33] A. V. K otikov and L. N. Lipato v, Nucl. Phys. B 661, 19 (2003) [Erratum-ibid. B 685, 405 
(2004)] [arXiv:hep-ph/0208220| . 

[34] A. V. Kotikov and L. N. Lipatov, [arXiv:hep-ph/0112346 

[35] A. V. Kotikov, L. N . Lipatov and V. N. Velizhanin, Phys. Lett. B 557, 114 (2003) 
|arXiv:hep-ph/030102l] . 

[36] A. V. Kotikov, L. N. Lipatov, A. I. Onishchen ko and V. N. Velizhanin , Phys. Lett. B 595, 
521 (2004) [Erratum-ibid. B 632, 754 (2006)] [arXiv:hep-th/0404092] . 

[37] A. V. Kotikov and L. N. Lipatov, Nucl. Phys. B 769, 217 (2007) |arXiv:hep-th/0611204] . 

[38] A. V. Kotikov, L. N. Lipatov, A. Rej, M. Staudacher and V. N. Velizhanin, J. Stat. Mech. 
0710, P10003 (2007) [arXiv:0704.3586ll [hep-th]]. 

[39] T. Lukowski, A. Rej and V. N. Velizhanin, Nucl. Phys. B 831, 105 (2010) |arXiv:0912.T624l 
[hep-th]]. 

[40] V. N. Velizhanin. larXiv: 1003.47171 [hep-th], 

[41] I. Balitsky and G. A. Chirilli, Phys. Lett. B 687, 204 (2010) [ arXiv:0911.5i"92l [hep-ph]]. 
[42] I. Balitsky and G. A. Chirilli. larX~iv:1009.4729l [hep-ph]. 



33 



